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ABSTRACT 

We address the issue of relativistic stellar aberration requirements for the 
Space Interferometry Mission (SIM). Motivated by the importance of this issue 
for SIM, we have considered a problem of relative astrometric observations of two 
stars separated by angle 9 on the sky with a single baseline interferometer. While 
a definitive answer on the stellar aberration issue may be obtained only in nu- 
merical simulations based on the accurate astrometric model of the instrument, 
one could still derive realistic conclusions by accounting for the main expected 
properties of SIM. In particular, we have analysied how the expected astromet- 
ric accuracy of determination of positions, parallaxes and proper motions will 
constrain the accuracy of the spaceraft navigation. We estimated the astromet- 
ric errors introduced by imperfect metrology (variations of the calibration term 
across the tile of interest), errors in the baseline length estimations, and those 
due to orbital motion of the spacecraft. We also estimate requirements on the 
data sampling rate necessary to apply on-board in order to correct for the stellar 
aberration. We have shown that the worst case observation scenario is realized 
for the motion of the spacecraft in the direction perpendicular to the tile. This 
case of motion will provide the most stringent requirement on the accuracy of 
knowledge of the velocity's magnitude. We discuss the implication of the results 
obtained for the future mission analysis. 

Subject headings: astrometry; techniques: interferometric, SIM; methods: ana- 
lytical; solar system; relativity 



Introduction 

One of the best promising methods of searching for a planetary systems around nearby 
stars is the high accuracy astrometric observations. SIM is being designed to produce wealth 
of the astrometric data necessary to address exactly this problem. It is clear that the 
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corresponding astrometric signal due to a reflex motion of a target star will be at the level of 
a few /ias and smaller. The most successful searching techniques will have to incorporate an 
intelligent way of data processing almost at the astrometric noise level. Another important 
part of this puzzle is development of the astrometric model for the instrument. This model 
should include all the parameters necessary to account for different phenomena affecting the 
light propagation, and are due to the interstellar media, the solar system dynamics, as well 
as due to the motion of the free-flying interferometer itself. As a result, the accuracy of 
astrometric observations expected with SIM, will require a number of dynamical parameters 
to be precisely known. To do this, one will have to be able to remove (or to correctly account 
for) the signatures of all the known effects in order to study the unknown phenomena. This 
is necessary for minimizing probability of a false detections caused, for example, by aliasing 
the dynamics of objects in the solar system. Actually, there will be a number of parameters 
introduced by the dynamics in the solar system with periods of the order of a few years. One 
of such a parameters is the three-dimensional vector of baryccntric velocity of the spacecraft. 
Thus, a future astrometric model will have to account not only for the effects due to the 
motion of the solar system bodies, but also for those that are generated by the motion of 
the spacecraft and corresponding errors in the spacecraft's navigation. 

In this paper we will discuss some important elements of dynamical model for the future 
astrometric observations with SIM. Speciflcally, we will address the issue of the relativistic 
stellar aberration. The stellar aberration is a very important problem for SIM. Effect of 
the barycentric velocity of the spacecraft is the largest term in the SIM astrometric model 
and will amounts to ~ 20.5 arcsec for the Earth-trailing orbit. This effect is important not 
only for a wide angle astrometric observations, but it will produce a measurable effect for a 
pairs of widely separated stars even inside the tile (with diameter ~ 15°). Additionally, a 
possible correlation of the astrometric obscrvables with the errors in the spacecraft's velocity 
sky-angles suggests that in order to achieve the mission accuracy for a global astrometric 
observations of cr = 4 /xas, one will have to control the barycentric velocity of the spacecraft 
throughout the entire mission and account for the relativistic stellar aberration inside every 
single tile. 

The outline of this paper is as follows. In Section 1 we will present a simplified model 
for differential astrometric observations with a free-flying interferometer with a single base- 
line. We will discuss a model for absolute and differential astrometric measurements and 
will give a three-dimensional parameterization of main observables. We will present the flrst 
order differentials, necessary to analyze the errors propagation. In Section 2 we will discuss 
orientation of the tile in the SIM nominal observing direction. We will present the first order 
equations that are governing the propagation of astrometric errors in the instrument. Specif- 
ically, we address the velocity knowledge requirement, the data sampling rate and orbital 
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position knowledge, and the spacecraft acceleration accuracy constraint. In Section 3 we will 
derive the set of the relativistic stellar aberration requirements for SIM. We will conclude by 
presenting our recommendations for future attitude control for the SIM spacecraft. In order 
to make access of the basic results of this paper easier, we will present some important (but 
lengthy!) calculations in the Appendices. 



1. A model for the relativistic stellcir aberration 

In this Section wc will derive a model necessary to analyze the astrometric error budget 
for differential astrometric observations with a single-baseline interferometer. To do this, we 
will present estimates of the optical path difference (OPD) for a singlc-basclinc interferometer 
in solar orbit. Our derivations will be different from the ones obtained earlier by the fact that 
for each tile we account for effects of general three-dimensional motion of the interferometer, 

— * 

for the errors in the components of the baseline vector b, and for the error in estimating 
the instrumental offset (or, calibration) term, cq. Our model is very simple, but it could be 
easily expanded to accommodate a number of other important features of the instrument. 

We address the problem of relativistic stellar aberration from a general standpoint. A 
"toy" astrometric model developed here will include only the largest relativistic contributions 
due to the orbital motion of the interferometer around the sun. A complete, fuUy-relativistic 
model for the SIM observations, will have to account for a number of physical phenomena 
and should include a set of additional terms (and corresponding parameters) necessary for 
different astrometric applications. More specifically, the future model should include the 
terms oc (^)^ to account for higher orders of relativistic stellar aberration; terms oc GAi/c^ 
(where G is the gravity constant and Ai is the mass of the deflecting body) to account for 
the gravitational deflection by the bodies of the solar system (primarily the Sun and the 
Jupiter). (For general discussion of the general relativistic effects in the SIM asdtrometric 
observations, see Turyshev (2002)) One also will need to account for rigid-body rotational 
motion of the spacecraft, etc. All these effects, together with a number of others, are out of 
the scope of the present paper, but will be discussed elsewhere. While the derivation of the 
general expressions will be given in the reminder of the document, here we will present only 
the main results obtained. 
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1.1. Geometry of the problem 

Let us begin by deriving the expression for OPD i for an interferometer with a single 
basehne h which is moving with the velocity f <^ c (for SIM this ratio will of of the order f/c ~ 
10~^) with respect to some reference frame (RF), for example the solar-system barycentric 
(SSB) reference frame. In the interferometer's proper quasi-inertial reference frame OPD, 
may be given as follows: 

e = {9 ■ s') + c'„ (1) 

— * — * 
where b' is the baseline vector, s' is the direction to the observed source, and the last term 

in this expression is the instrument offset (or "constant" term), cq, that can be calibrated 

out. The same quantity may be presented in the coordinates of the SSB reference frame. 

The corresponding transformation between the solar system barycentric frame and the one 

of the moving interferometer (with velocity v with respect to SSB) will necessarily account 

for the effect of the relativistic stellar aberration (for more details see Turyshev (1998)). To 

the first order in ^, the aberrated direction to the source and the baseline vector may be 

expressed as follows: 

- ^ ,^-:_!(^ + 0(G,(^r), (2) 

c c V c / 

9 = b + o{G,{y), (3) 



where s and b are the direction to the observed source and the baseline vector as measured 
in SSB. The total effect of this transformation has the form 

e = (i.s,{i-^l^) + <t^ + c + o{G,C-r), (4) 

where we have assumed that the instrument offset does not depend on the velocity of the 
spacecraft motion to the second order in ^, or, in more general terms, Cq = cq + o(^G, (^)^^ • 
This issue may be correctly addressed later when the accurate model of the instrument will 
be developed. 

The simplified two-dimensional geometry of the problem presented in Figure 1. This 
allows to express the scalar products in Eq.(4) and re- write this equation in a more familiar 
form 

/ V \ bv 
£ — bcosail cos(q; — -0) I H cos -0 + 00. (5) 
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Fig. 1. — Two-dimensional geometry of the problem: absolute astrometry. 

Such an expression may be used to model absolute astrometric observations with a moving 
interferometer, when the observed sources are not in the same field of view, but rather widely 
separated from each other. However, the wide-angle astrometric campaign with SIM will be 
based on a set of differential (or relative) astrometric observations made inside a set of 637 
tiles covering the whole sky (Boden (1997); Swartz (1998)). As a result, it will measure to 
a certain accuracy the angular distances between a pairs of stars in the same field of view. 
This necessitates the derivation of an expression for the optical path difference, similar to 
Eq.(5), which would take this fact into account. To the first order in | this is a fairly easy 
task. Assuming that absolute position for each star is given by expression, similar to that of 
Eq.(4), namely: 



one could write the corresponding differential (or relative) OPD for a differential astrometric 
observations. To describe differential measurements one will have to subtract one OPD from 
the other, namely 81 = l\ — l2- The corresponding differential OPD will be given as 
follows: 




(6) 




S\ ■ V 



c 



S2 ■ V 



c 



) 



(7) 



where 6co is the differential "constant" term 6co = cqi — co2- 
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1.2. Three-dimensional parameterization. 

In this subsection we wiU present the three-dimensional parameterization of all the 
vectors involved in the problem. This from of representation is more useful for practical 
applications, for example, to study the orbital motion of the spacecraft; for analyzing the 
various requirements on the spacecraft's velocity, rigid-body rotations, vibrations of the 
whole structure, etc. 

The overall three-dimensional geometry of the problem presented in Figure 2. The 
different vectors involved are given in the spherical coordinate system by their magnitudes 
and two corresponding sky-angles: 

si = (cosai cos(5i, sinaicos^i, sin^i), 

S2 = (cosa2Cos52, sin 02 cos (^2, sin 52), 

V = f ( COS t/' COS e, sin '?/' cose, sine), 

b = 6( cos/x cos z/, sin /i cos z/, sini/). (8) 



- 7- 



This parameterization allows one to present expression Eq.(7) in the following form: 

sin di — sin 62^ 



5i = b 



bv 
c 



C0Si/^C0s5iC0s(Q;i — //) — COS(52COs(q;2 , ,j 



sin 61^ ^ cos e cos 5i cos(q;i — tjj) + 



cos u cos ^1 cos(q;i — /x) + sin u 
^ cos 1/ cos S2 cos(q;2 — /i) + sin u sin ^2 j ^ cos e cos ^2 cos(q;2 — ■0) + sin e sin ^2 



+ Sco - 
sin e sin 5i ) — 



(9) 



This expression is quite difficult for analytical description, however, for the purposes of 
the present study, it may be simplified (a more detailed analysis of this problem will be given 
in the Appendix A). Thus, one may see that the angles of the baseline orientation (/x, i/) 
are significantly influencing the narrow-angle astrometric observations. Angle fi corresponds 
to a shift of the origin of RA for the tile under consideration and, without lost of generality, 
it may be omitted (e.q. set to be zero). The analysis of contribution of the DEC angle u 
is a little bit more complicated task. For the purposes of this paper we will assume that 
this angle will he i/ — at the beginning of the tile observation and will insignificantly 
deviate from this value during the experiment. This is true, because one of the functions 
for the two guide interferometers in SIM will be to provide a stable reference for the science 
interferometer. We, therefore will set both angles as ^ = z/ = 0. These parameters will 
have to be incorporated in the future astrometric model for SIM, probably in a form of 
the rigid body rotation/precession/nutation model of the spacecraft. At this moment, this 
choice is equivalent to choosing the direction of the baseline vector b coinciding with x-axis: 
b = 6(1,0,0) (see Figure 2). As a result, all vectors now will be counted from the baseline 
and the expression (9) may now be presented in a simpler form, namely: 



b cos 61 cos ai — cos 62 cos q;2^ + Scq — 
^ ^ cos 5i cos Oil ^ cos e cos 5i cos(q;i — -0) + sin e sin 
■ cos 82 COS a2 ^ cos e cos 82 cos(q;2 — ■0) + sin e sin ^2^ 



(10) 



The obtained expression may now be used to analyze the error propagation for the 
moving interferometer. To do this, one will have to expand this expression in a Taylor series 
and keep only the first term of the expansion. While the details of this calculations are 
given in Appendix A, here we present only the final expression which is necessary to analyze 
contributions of the chosen set of different error factors to the overall error budget: 
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A5£ 



= Aq! sin Q!2 COS 62 + Aq!i ^ sin a2 cos ^2 — sin ai cos + 

+ A5 cos a2 sin 5^ + A5i ^ cos a2 sin ^2 — cos ai sin + 
Ab f , N A5co 

+ — — cos Oil COS Oi — COS a2 cos ()2 I H ; h 

V / b 



+ 



— cos(e — eo) — Ae - sin(e — eo) \/ (a^ + p) cos^lip — ipo) + — 
c c J 



Alp - cos ei/ 0? + /2 sinf'^ — ■i/^o)) 
c 



where we have used the following definitions for the situation discussed: 



a — a2 — Oil 
5 — 82 — 5i 



Aa 
AS 
Aai 



A5i 



Ab 



AS, 



Co 



Av, Ae, Ai/j 



RA angle between the two stars inside the tile under consideration; 

DEC angle between the two stars inside the tile under consideration. 

[Note that the following identity is expected to hold for these two 

angles + 5^ < {^f ]; 

error in defining a: Aa = tttrue - ^estimated; 

error in defining 5: A5 ^ Stme - (^estimated; 

error in defining RA angle, ai, for the primary star 

^^^1 '-^Itrue '-^l estimated 7 

error in defining DEC angle, 5i, for the primary star 

A5i — ^Itrue ~ "^1 estimated) 

error in defining the baseline length b: Ab = b-^rue — ^estimated! 
error in defining the relative constant term for the observations 
inside the tile: AS, 



CO 



■'CO true 



CO estimated' 



errors in knowledge of components of the three-dimensional spacecraft 
velocity vector, defined as usual =^ A = true - estimated. 



Additionally, the new quantities ipQ and are entirely defined by the coordinates of the 
two stars under consideration and are given as follows: 
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, , COS'^o = , , 

COS^ ^2 COS^ (X2 — COS^ COS^ OL\, 

cos^ ^2 COS q;2 sin ol^ — cos^ 5i cos a\ sin cci. (12) 
k 

■y/g^ + cos(V> - V>o) 
v/(a2 + /2)cos2(V'-Vo) + A;2' 

COS ^2 sin ^2 COS 0:2 — cos 8\ sin 5i cos cti. (13) 

These expressions (11)-(13) may now be used to study propagation of astrometric errors 
to the first order. With these results one may analyze the relativistic stellar aberration issue 
analytically for a different positions of the two stars of interest. Note that in reality all the 
quantities involved are both dynamical and stochastic functions of time. This fact we be 
used in the Section 3 when we will analyze the allocations adopted for different constituents 
of the SIM astrometric error budget. 

2. Tile in the SIM nominal observing direction 

Analysis of the tolerable errors for the velocity aberration in a general case is a compli- 
cated problem and should be attacked by using a formal numerical treatment. However, we 
may simplify the task by analyzing a number of a special cases which are expected to provide 
the most stringent requirements on the velocity magnitude. Thus, one may expect that the 
most driving requirement will come in the case when the two stars are aligned in a tile in a 
such a way that the angular separations between them is given by a = 5 = 0. It turns 
out that this is indeed the case. The compensation for aberration across the field of view of 
the interferometer in a direction parallel to the baseline, will present the greatest challenge 
for SIM. To demonstrate this, we will discuss a more general situation (that includes the 
case mentioned above), which allows one to obtain an analytical expressions useful for the 
future analysis. 

Let us define a tile in the SIM nominal observing direction — perpendicular to the 
baseline. For the purposes of our analysis we will assume that positions of the two stars 
are symmetric with respect to the point with coordinates (cto = f ; = 0). Then the most 



sin — 

a — 

f = 

sin eo = 

cos eo = 

k = 
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Fig. 3. — A tile in the SIM nominal observing direction. Note that + 5^ < (^)^. 



suitable description for the nominal SIM observing mode is given by coordinates of the 
primary and the secondary stars as follows: 

71 a TT a 

2~r "^=2 + 2' 

5i = 52 = ^. (14) 

Now we can define a tile for this observing direction as a set of all the points for which the 
angular separations a and 6 are limited as: a'^ + S'^ < (f^)^. The resulted area constitutes the 
tile in the SIM nominal observing direction and it is shown as a shaded area in the Figure 3. 

By substituting the coordinates for the two stars Eq.(14) in the expressions (12)-(13) 
we obtain the following values for the constants involved: 



a 

f 
k 



sin a cos^ |, 



• 2 4 <5 

sm a cos 
sin^ a cos^ |. 



sin i/jq 

COS'?/'o 

sineo 
coseo 



-1 



1 



2 ' 



eo = 0. 



Utilizing these results in the equation (11), we obtain the expression necessary to analyze 
contributions of different error factors to the total differential OPD. This expression has the 
following form: 



A5£ ^ a 5 
—— — Aa cos — cos - 
h 2 2 



... a . 5 2Ab . a 5 AScq 
Ad sm — sm - + — — sm — cos - + — ; — 
2 2 6 2 2 b 



fAv ^ ■ \ ■ I X I ^ 

— cos e — Ae - sm e ) sm + A'^ - cos e cos ip 
\ c c J c 



sm a cos 



(15) 



The obtained expression Eq.(15) may now be used to analyze the propagation of errors 
in the future astrometric observations with SIM. 



2.1. Propagation of cistrometric errors 

In this section we will concentrate on obtaining the expressions that are essential for 
derivations of the requirements on the quality of the spacecraft navigation data. These 
requirements are imposed by the expected accuracy of the future SIM astrometric measure- 
ments, such as determination of positions, parallaxes and proper motions. 

In order to derive the necessary equations, the first term in the expression (15) A5£ 
may equivalently be presented as A6i = A(nAo) = AnAo + n AAq, where Aq is the operating 
frequency and n is an integer number. This term vanishes because both Aq and n are assumed 
to be known to a sufficiently high accuracy, thus An = 0, AAq = 0. This is true because 
SIM will be using it's fringe tracker to find a white light fringe to perform the calibration of 
the instrument. Then the remaining differentials Av, Ae, A^, Aa, A6, Ah and AScq will 
have to satisfy the equation: 

. a S A r A6co 2Ab a 5 

Aa cos — cos - — Ad sm — sm - = ; ; — sm — cos — h 

22 22 6622 



+ 



— cose — Ae - sme sm'^ + A'^ - cose cos smacos -. (16) 
. c c / c J 2 



S 



In order to simplify further analysis wc have separated the terms in the expression above in 
a such a way, so that the left side of this equation represents the error in the measurement 
of absolute angular separation between the two stars on the sky, while the left side shows 
the main contributing factors to this quantity. 
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2.1.1. Velocity knowledge requirement 

One may expect that, in any given tile, the errors Aa and A6 are normally distributed 
and uncorrelated.^ Also, the errors due to the orbital dynamics Av, Ae, Aip at the chosen 
approximation may be treated as being uncorrelated with the instrumental errors Ab, Acq.^ 
Then, one obtains: 

2 2 « 2 <^ , 2 • 2 « • 2 ^lo , . 2 a 2 ^ , 

(7„ COS - COS - + as sm - sm - = + ^ 2 

/ O^q, 9 2^ •2^*2/ 2^ 2 2/1-2 4^ / 1 

+ [—^ cos e + (7g sm e I sm + cr^ — cos e cos sm a cos -. (17) 

For the wide-angle astrometric observations with FoR~ the interferometer is ex- 
pected to perform at the level of cTq, = 4 fias and ci^ ~ 1 mas. Expression (17) suggests 
that for observations in the direction perpendicular to the baseline (e.q. when angle 6 takes 
non-zero values), there will be a large contribution coming from as to the error budget, 
thus weakening the relativistic stellar aberration requirement. However, for the observations 
parallel to the basehne (e.q S — 0), influence of as is vanishes and one obtains the following 
equation: 

2 2 ^5cn 4cr? . 2 Oi 

''^'°' 2 = # + 1^™ 2 + 



+ 



^2 „,2 , „,2 



{^-^ cos^ e + a'^ ^ sin^ e j sin^ + a'^ — cos^ e cos^ -0 sin^ a. (18) 



This expression depends on the orientation of the spacecraft velocity with respect to 
the SSB coordinate reference frame. Let us analyze the worst case of the spacecraft orbital 
orientation. This orientation is when the spacecraft is moving towards/from the tile in the 
direction exactly normal to the tile of interest. In this case, when ip = ±|, e = (or in the 
direction normal to the tile), from the equation Eq.(18) one immediately obtains expression 
that contains a requirement on the error in the magnitude of the spacecraft's velocity: 

2 na cr? 4(j2 Q, ^2 
'""'^^' 2 ^ # + "6^ '^" 2 ^^^^ 



^Note that this is not true for a general case of studying the stability of the reference grid (Swartz (1998)). 
Thus, one finds that the correlation in RA and DEC becomes a source for the zonal errors in the analysis of 

the grid accuracy and stability. 

more general analysis should include a possible correlation between the constant term Acq and the 
two angular components of velocity Ae, A'^. We will address this possibility later in the Appendix C and 
will discuss the components cov(co, e) and cov(co, V') of the total covariance matrix. 
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One could also verify, that errors in both sky angles t/j and e are related to that for the 
velocity magnitude v and are given by 



^2 ^2 "^v 



(20) 



Note that a possible correlation between the constant term 5cq and the two sky angles 
of velocity ip^ e seriously affecting these estimates (a more detailed analysis of this problem 
is given in Appendix C). These errors, and a^, are related to the error in the velocity 
magnitude as follows: 



f O'^co P(C0,6) CI 2 

b sina v- 



+ 



^^CQ p(co,V') c 
h sina v 



o"^c„ p(co,e) c 
6 sin a V 

6 sin q; i> 



> 0, 



> 0, 



(21) 



where 



(7e, (7^ - the errors in the two sky angles of the bary centric velocity vector v] 
p{co, e), p(co, ■0) ~ correlation factors between the constant term and the two sky angles 
of the spacecraft's barycentric velocity, |p(co,e)|, \p{cQ,ip)\ < 1. 

In addition to the velocity knowledge requirement one will have to impose two additional 
constraints on the quality of navigation data delivered by DSN. These two constraints are on 
the sampling rate At and positional accuracy of the spacecraft determination Ar together 
with the requirement on the stochastic acceleration control during the time of astrometric 
observations. 



2.1.2. Sampling rate and orbital position knowledge 

As we discussed previously, the correction for the annual stellar aberration may reach 
~ 20.5 arcsec. This allows estimation of the time interval on which this correction must 
be introduced in order to maintain the nominal accuracy. Thus, contribution of only the 
velocity of the spacecraft to the total error budget Aa may be obtained from Eqs.(16). For 
the motion in the plane of ecliptic with zero declination (5 = e = 0), one will have: 

. a /Av . , J, , V A . 

Aa cos — = — sm'^ + Aip - cos'^ sma. (22) 
2 V c c / 



How often one will needs to introduce a correction for the stellar aberration? This 
question may be answered directly by analyzing the motion of the spacecraft in the plane 
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parallel to the tile, when the aberration takes it's largest value. The latitude argument 
for circular motion coincides with the mean anomaly, e.q. i/j = nt, where n = ^ and 
P is the period of orbital motion of the spacecraft. Taking this into account and setting 
(-0 = 0; e = 0), from equation (22) one obtains: 

/dAa\ a ^ , v . 2-k ^ , „, 

— — cos- Alb = - sm a — At, (23) 
\ dip J 2 ^ c P ' ^ ' 

where A'^ — nAt and At is the time interval. After averaging, the last equation may be 
presented in a more useful form (note, (Tq ~ Atp), namely: 

(Tq, cos — > - sin a -5- At, (24) 



As a result, one will have to correct for the relativistic stellar aberration with such a sampling 
rate At that the expression on the right-hand side of the equation above will be considerably 
less then on the left side. 

The data samphng rate At may be used to define the corresponding accuracy of de- 
termination of the spacecraft's position on it's orbit Ary. By decomposing the spacecraft's 
orbital position vector on the component parallel to velocity vector and the one perpendicu- 
lar as given by Af— Af± + Afji, one may estimate the magnitudes for the both components 
involved. The expression for the magnitude of Afji may be obtained simply as follows: 

Ar\\<VsmAt. (25) 



A requirement on the magnitude of the second component of the spacecraft position 
vector Ar_L may be derived from another mission goal, namely the expected accuracy for 
parallax determination. Thus, SIM is expected to achieve a mission accuracy of the de- 
termination of parallax at the level of cr^ = 1 /xas for a nearby stars. This goal may be 
transformed into the requirement on the accuracy of the radial component of the barycentric 
spacecraft position Afj_. On may show that the error, introduced by annual parallax into 
the astrometric observations of a pair of stars separated on the sky by angle a (with S = 0), 
may be given as below: 

r r^l\ ( I oii + a2\ . 02 -ai 

6a^ai-a2 = —2 cos [ip j sm , (26) 

where r|^' is the spacecraft barycentric distance (the superscript [tt] is used to separate 
parallactic distance requirement from those derived with the help of some other methods), 
and V is the distance to the object of study. Then for the tile in the SIM nominal observing 
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direction, defined by Eqs.(14), and for the motion perpendicular to the tile {ip — ±|) one 
may obtain following relation: 

Ar!r^< a.^^-^. (27) 
2sm| 

This relation suggests that the closer the distance T> to the observed object the better should 
be the knowledge of the barycentric position of the spacecraft. 

Accuracy of the spacecraft orbital position determination Ar is one of the navigational 
products that will be provided by means of DSN. However, one should not expect that DSN 
will provide such a continuous trucking of the spacecraft with the data sampling rate At. This 
constraint (primarily cost of such an extended DSN commitment) implies that a real-time on- 
board correction must be done in order to provide the conditions necessary for astrometric 
observations with expected accuracy. Note that this problem may be already solved by 
means of the SIM multiple-baseline architecture, when the two guide interferometers will 
take care of a number of similar effects. In any case, this question is out of scope of the 
present analysis, but it should be re- visited as a subject for separate study. 



2.1.3. Acceleration knowledge constraint 

The two additional constraints, namely on the quality of the DSN data indirectly im- 
posed by the expected accuracy for proper motions determined for a number of stars during 
the mission's five-years live time, which will be of the order of /Xq, = 1 /xas/yr. Such a con- 
straint may be also obtained from the Eq.(22). By differentiating this equation with respect 
to time one obtains: 

dAa a /I dAv . , dAib v 
cos ■ ' 



sm'^ -|- — ; cos-^lsmQ;. (28) 

dt c / 



dt 2 \c dt dt 
The time derivatives on the right-hand side of this equations are the components of the 
spacecraft acceleration. Indeed, one may decompose vector of this acceleration onto the 
components normal to the orbit (radial) and the tangential one as follows a — d± + d\\, 
where = a-nd ay = Assuming that the errors in both accelerations are normally 
distributed and uncorrelated, this equation may be averaged and presented in the following 
useful form: 

2 7,2 



2 2 
a,, cos — 

Ma 2 



For the motion of the spacecraft in the direction normal to the tile of interest (V' = ±f ), 
one may obtain expression that imposes a requirements on the knowledge of the spacecraft 
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radial acceleration the following form: 



a 1 (iAf . , . 

<^Ma cos - > — ^ sm a. (30) 
I c at 

Similarly, for the motion of the spacecraft in the direction parallel to the tile {ijj = 0), one 
will have a requirement on the knowledge of the angular acceleration: 



a V dAih 

(J^a cos- > --^sma. (31) 



These two equations suggest that, for accurate determination of proper motions one 
will need to control all the non-gravitational forces acting on the spacecraft at a very high 
level of accuracy. Note, for the expected orbital parameters of SIM, the control of the radial 
accelerations should be at the level of ~ 1.2 x 10~^^ km/s^. 



3. Relativistic stellar aberration requirements 

In this Section we will present the estimates for a different strategies of accessing the 
influence of the relativistic stellar aberration on the expected accuracy of astrometric mea- 
surements with SIM. 

We will begin from the expressions describing the uncorrelated data Eqs.(19)-(20). Re- 
membering that angles a, S vary in the range given by + < (f^)^ for which cos | never 
vanishes, one can divide the both sides of the equation (19) on cos 



2 a. 
2 ■ 



The obtained equation represents ellipsoid with half-axis defined as follows:^ 

o-5eo ^ f^a&cos|, (33) 
h a 

o"6 = ^a2Cot-, (34) 
2sm I 



more general expression for the stellar aberration only may be obtained from Eqs.(ll)-(13). It depends 

.(ct2— « 

sin 02 



on the coordinates of both stars under consideration and has the following from: Gq^^ = — ^"^("^ "i) 
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The necessary expressions for the tolerable errors in the components of the three-dimensional 
vector of the spacecraft's velocity may be obtained directly from Eqs.(20),(35) as 



^^^^ 

In addition to these four constraints we will have to define the other three discussed in 
the previous Section 2.1, namely 1) the data sampling rate, At; 2) the accuracy of the orbital 
position determination, Ar± and Ar||; and 3) the accuracy of compensation of the non- 
gravitational accelerations, a± = and a\\ = By utilizing equations (24)-(27) together 
with Eqs.(30),(31), these constraints may be presented by the following set of equations: 



dAv c dAi) c 1 

< a,, , < a,, . (39) 

dt - '^"2sinf' dt - '^"vsi„2sin" ^ ^ 



These five requirements arc the derived constraints on the orbital motion of the spacecraft 
and the data sampling rate. While the stellar aberration constraint is the most important 
among those imposed by the orbital dynamics, we think that it is useful to have all the 
requirements related to the motion of the spacecraft spelled out. 

Let us assume that a different error-budget constituents are contributing individually 
to the total error budget with an individual astrometric error cr^ = Ak /xas, where Ak is 
some number. By taking into account the numeric values for different quantities involved 

(7« = AA; //as, c = 2.997292 x 10^^ mm/s, 

b = 10.50 m, 1 //as = 4.84814 x 10~^^ rad, 

(7^„ = Al //as/yr, a^, = Aq //as, 

= 2.98 X 10^ mm/s, P = 3.1536 x 10^ s (40) 



Vsm 



one could compute tolerable astrometric contributions introduced by the different errors 
budget constituents. The largest effect (the worst case) will be in the case of maximal 
separation between the two stars which will equal to the field of regard of the instrument 
Qimax =FoR= j^. Substituting the numerical quantities from (40) we immediately obtain the 



following values for various components of the error budget: 



as^^ = 50.470 A/c pm, 
(76 = 193.332 AA; pm, 



(7v — 5.566 AA; mm/s, 



— ae — 38.529 Ak mas; 



(41) 



An 



II 



27.952 AA; km, 



At = 0.938 AA; s, 

Arf ' = 573 Aq km. 




dAv 



11.132 Al mm/s/yr. 



(42) 



dt 



where we have assumed that the barycentric velocity of the spacecraft is approximately equal 
to that of Earth' orbital motion around the sun, (e.g. c/vsm = 10058), thus giving us the 
requirement for the two sky angles of the velocity vector. 

The constraints on the spacecraft's position are given by the Eqs.(42). It is worth noting 
that the parallactic requirement on the spacecraft barycentric position Ar|[^' is easy to meet 
for a distant objects, because it grows linearly with distance. [Note that for the observa- 
tions of the solar system objects one will have to use a completely different observational 
strategy.] The accuracies of the DSN navigation are much superior than it needed to sat- 
isfy this parallactic requirement. Note that the general relativistic deflection of light is also 
distance-dependent effect. One may think that this dependence may produce an independent 
constraint on Ar||^'. However, a crude estimate shows that this is not the case and general 
relativity does not require a significant accuracy of knowledge of the barycentric distance of 
the spacecraft. This is why we will omit the constraint on Ar^^^ from our studies. The last 
two constraints in Eqs.(42), on the two components of the spacecraft acceleration. These con- 
straints are suggesting that the total uncompensated error in the acceleration of the craft over 
a half-orbit (half-year time interval) of the mission should not exceed these numbers. The 
first requirement may be re- written in more familiar units as = 3.53 x 10~^^ Al km/s^, 
which again implies quite a significant navigational involvement or a real-time on-board 
processing. 

The worst case observation scenario will be realized for those angles ip and e which will 
maximize the corresponding multipliers in formulae (18). Thus, the motion of the spacecraft 
in the direction perpendicular to the tile, or {ip ^ ip^ = e — 0), will provide the 
most stringent requirement for the accuracy of knowledge of the magnitude of barycentric 
velocity of the spacecraft. Additionally, the motion in the direction parallel to the tile. 



-19- 



to 



Pi 



\ / 

/ / 




\A / \ 





Spatial fringe resolution drift 

Fig. 4. — Spatial fringe drift introduced by the relativistic stellar aberration. Suppose the 
data are taken with the sampling rate At, then the thick curve - is the fringe intensity at 
the beginning of observations at time Iq, while the dotted curve is plotted for the position 
of the same fringe after time interval /St has elapsed. 



or {ip — > ipw = 0,7r; e = 0), will provide the most stringent requirement for the accuracy 
of knowledge of the angular position of the spacecraft with respect to the solar system 
barycenter and, therefore, it introduces constraint on the data samphng rate At. This rate 
is necessary to apply for on-board navigation in order to correct for the stellar aberration. 
The nature of this correction for a single-baseline interferometer is presented by the Figure 
4. 

Relativistic stellar aberration introduces a spatial fringe drift due to the motion of the 
craft. This spatial fringe drift may be corrected if one will be sampling data with the rate 
At. This time interval tells how often one will need to introduce a correction to the velocity 
of a single-baseline interferometer in order to provide the necessary dynamical conditions for 
astrometric measurements with accuracy (Jq, (i.e. for the starts separated by FoR=«). SIM 
is the multiple baseline interferometer. The main function of the two guide interferometers 
is to provide this sort of anchoring to the sky by providing a set of differential navigation 
parameters. Moreover, for the direction of motion parallel to the tile of interest, the accuracy 
of the angular position of the spacecraft is more important then it's velocity magnitude, thus 
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relaxing quite a bit the corresponding values for Ak in Eqs.(42). So, to the first order, the 
requirement imposed by At does not significantly impacting SIM astrometric campaign, but 
it is important enough to consider it for further studies. 

3.1. Two ways of direct establishing the error allocations 

In this Section we will present analysis of different methods for establishing the errors 
allocations for the stellar aberration. The expressions Eq.(41) are representing the set of 
requirements (including the stellar aberration) for the astrometric observations with accuracy 
of (Tq, = Ak /ias in the direction parallel to baseline. By varying the coefficient Ak wc will 
generate the sets of requirements for as^^, Cb, ^V' '^e \hsX correspond to these different 
methods of error allocations. 

There exist a mimbcr of different ways how to estimate contributions of different con- 
stituents to the total astrometric error budget. For a more accurate analysis, one should 
perform a numerical simulations in order to estimate the impact of those constituents on the 
astrometric grid performance. At this point we would like to analyze a number of possible 
limiting cases, in order to provide a taste of different flavors to this problem. To do this 
we will use the results presented by the set Eqs.(41). In particular, we will analyze the 
constraints that should be placed on different error budget constituents in the cases when 
one choose to make the assumptions based on either single measurement accuracy or the 
expected wide-angle mission accuracy. 

3.1.1. Single observation accuracy requirements 

SIM is being designed to be able to achieve astrometric accuracy of a single measure- 
ments of the order of (Tq = 8 //as. This is first useful number to use, in order to derive the 
maximal stellar aberration error, tolerable during the tile observation. The corresponding 
estimate may be obtained by substituting the value Ak — 8 into Eqs.(33)-(35). As a result. 
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we obtain the following values for different terms in the equation Eq.(32): 



,[81 



[8] 



403.76 pm, 

1546.66 pm, 

44.53 mm/s, 

308.23 mas; 



(43) 



dAv 



[/\l= 



= 44.53 mm/s/yr, 



7.50 s, 
223.62 km, 

;[A/=8] 

= 616.46 mas/yr 



dt 



(44) 



A reasonable improvement, for the numbers presented above, comes from the statement 
that a contribution of any component of the total error budget in the right-hand side of the 
equation (32) should not exceed 10% of the total variance a single accuracy of o"^. This gives 
a different correction factor AA:°'^ which is calculated to be AA;°'^ = Sv^OT = 2.52982. The 
resulting numbers for the contributions to the error budget in a single measurement mode 
are follows: 



a 



J2.5] 



[2.5] 
J2-5] 

V 

[2.5] 



dAv 
~df 



[Al=2.5] 



14.08 mm/s/yr. 



127.68 pm, 

489.10 pm, 
14.08 mm/s, 
97.47 mas; 

2.37 s, 
70.72 km, 

-[A/=2.5] 



(45) 



194.94 mas/yr (46) 



It is important to point out that the velocity aberration issue has two flavors to it, in 
a sense that it is influencing both the narrow angle observations and (by the procedure of 
the astrometric grid reduction) the wide angle ones. This mean that the error coming from 
the velocity aberration knowledge inside one single tile, will propagate to the total accuracy 
of the wide angle observations. This is why it is important to estimate the corresponding 
maximal tolerable errors based on the expected wide angle astrometric accuracy. 

The set of requirements represents Eqs.(45),(46) an optimistic expectation on the ac- 
curacy of the velocity determination. Due to the reason that the astrometric accuracy will 
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be improving as mission progresses, this set of relativistic stellar aberration requirements 
will be easily met by the DSN navigation. However, the assumption based on the single 
measurement accuracy is over-optimistic and we need to consider the most driving cases 
presented by the expected wide-angle mission accuracy 



3.1.2. Mission accuracy requirements 



The wide angle astrometric observations with SIM are expected to be with a mission 

accuracy of do = 4 /xas. Using this number, one may derive a different set of requirements 
for SIM, which will be exactly twice smaller than the ones given by Eq.(43). Additionally. 
by applying a conservative 10% approach (described above), one may derive another set of 
estimates. The wide angle astrometric mode together with the conservative approach, gives 
for Ak the number A/c"'^ = 4\/0.1 = 1.26491. This value of the correction factor results in 
the following contributions to the error budget: 



[1.3] 



^[1-3] 
Jl-3] 
r[l-3] 

V 

[1.3] 



dAv 
"df 



[Al=1.3] 



= 7.04 mm/s/yr. 



63.84 pm, 

244.55 pm, 

7.04 mm/s, 

48.74 mas; 

1.19 s, 

35.36 km, 

-[AZ=1.3] 



dAl/j' 

~dt 



— 97.47 mas/yr. 



(47) 

(48) 
(49) 
(50) 

(51) 
(52) 

(53) 



There is another important issue which hasn't been addressed yet. That is a possi- 
ble correlation between the calibration term 6c^y and the velocity sky- angles {ip, e) given by 
Eq.(21). To estimate the influence of this possible correlation, we assume that the two cor- 
relation coefficients for the two angles involved are equal p{co,ip) = p(co,e) = po- Then, 
by using the expressions (21) together with the estimates Eqs.(33)-(35) one could obtain an 
improved expression in order to derive requirements on the accuracy of knowledge of the two 
sky-angles of the spacecraft's velocity: 



a.^a^^aa^ (^^l + p^-po) = (^1 + - po) 



> 0. 



(54) 
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It is important to point out, that for the worst case of highly correlated quantities {ip, e) and 
Co, e.q. po 1 the result Eq.(50) will have to be further reduced by a factor of \/2 — 1 = 
0.4142, thus tightening the requirements for the accuracy of knowledge of the two velocity 
sky-angles as 0"^ = o"^ ~ 20 mas. [Note that this estimate is for the worst case of correlation, 
when p ~ 1. For the case, when po ~ 0.5 this correction factor is almost twice relaxed, namely 
(1 — 0.5^)^ — O.S'^ = 0.868.] This example suggests that a possible correlation between the 
velocity components and the constant term in the narrow-angle observations may put an 
additional demand on the quality of the DSN data. Thus for the worst case scenario, DSN 
will have to deliver data for all three components of the velocity vector with considerably 
smaller errors, say 

(7^°^^ = 2.916 mm/s, (55) 
cr"^^ = (7,=°^^ = 20.186 mas. (56) 



The two numbers above arc the very pessimistic estimates for the required accuracy of 
determination the spacecraft velocity vector and are given for the worst case of highly corre- 
lated data. We have obtained this requirement by using the number for the expected mission 
accuracy for wide-angle astrometric observations together with a conservative allocation^ for 
different errors in the total astrometric error budget. 



3.2. Currently adopted error allocations 

Presently the analysis of not only the stellar aberration, but also a number of the other 
important issues is complicated by the fact that a reahstic model for the spacecraft and the 
instrument is absent. This is why the currently adopted values in the error budget were 
chosen more or less intuitively. Moreover, these numbers are the same as they were for the 
Earth-orbiting mission study. Now, when SIM is designed to be placed on the Earth trailing 
SIRTF-like solar orbit, these numbers should be verified for a better justification. The reason 
for doing that is not only the need for justification of the number of mm/s (this number is 
the same for both orbits), but rather time and the level of DSN commitment, necessary to 
reach the desired accuracy of velocity determination. 

Remember that the numbers in Eqs.(50) were obtained based on the assumption that 
the errors on the right-hand side of the Eq.(19) are forming an ellipsoid with half-axes given 



"^This estimate was derived from the statement that the contribution of any component of the error budget 
in the right side of the equation (32) should not exceed 10% of the total variance cr^, with CTq = 4 /las. 
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by Eqs.(33)-(35). In reality, one will have to minimize each constituent of the total error 
budget in a such a way that in any given time the sum of the terms on the right-hand side 
of the equation (19) will not be larger then the expected variance a^. Currently, the error 
budget estimations allocates for the stellar aberration = 36 pm. The corresponding Ak 
is then estimated to be of order Ak = 0.701 [based on a^^ = (a^^ + f^sj^ ~ = = 36 
pm and 6=10.50 m]. The requirements now will have to be modified as follows: 

(57) 

(58) 
(59) 
(60) 

(61) 
(62) 

54.02 mas/yr (63) 



dAv 
If 



[AZ=0.7] 



Acq 


35.38 


pm. 


[0.7] 


135.53 


pm. 




3.90 


mm/i 




27.01 


mas; 




0.66 


s, 




19.59 


km, 


3.90 mm/s/yr. 


dAil) 
dt 


[Al=0.7] 



Despite the fact that we have a reasonable gap between our estimates Eq. (49) and the 
best experimental guess given by Eq.(59), there some other factors that are necessary to 
consider. As we demonstrated previously, a possible correlation between the constant term 
and the two velocity sky-angles Eqn.(54) may completely ehminate this gap and further 
reduce the estimates presented in the Section 3.1.2 (for the worst case of highly correlated 
data). This fact minimizes a tolerable errors Eqs.(47)-(53), reducing those down to the values 
Eqs.(57)-(63). 

Accounting for this correlation directly in the expressions above, leading to — 1.616 
mm/s. However, remember that the contribution of stellar aberration to the total error 
budget was chosen to be cr^^ = 36 pm. In order to see whether or not this number should 
be kept as a maximum tolerable stellar aberration error, this problem should be addressed 
with a formal numerical treatment. In the mean time, we recommend that a^^ be reduced 
to, say ~ 25 pm. 

A possible increase of the field of regard (FoR), which is currently being discussed, will 
also results in minimizing a^. In the Figure 5 wc present the dependence of the maximal 
tolerable error in the velocity determination function of FoR. One can see that 

relativistic stellar aberration increases roughly as the tile angle squared. Thus, increasing 
the tile size from 15 to 20 degrees puts a tighter requirement on the knowledge of the 
spacecraft velocity which is already difficult because of the earth traihng orbit. Therefore a 
further study of this problem is granted. 
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Fig. 5. — The size of the tolerable velocity error as a function of FoR. Plotted from Eq.(35) 
for the astrometric error of o"„ allocated for the relativistic stellar aberration in the total 
error budget. The upper dashed curve correspond to cTq, = 1 /xas and the lower thick curve 
is for (Ta = 0.701 lias. 



3.3. Requirements for on-bo£ird attitude control 

In this Section we will estimate the influence of a temporal changes introduced by 
the rotational motion of the interferometer with respect to it's center of mass. Thus, the 
expression (50) may be used to derive another requirement, namely the requirement on 
the knowledge of the rotational motion of the spacecraft, or on-board attitude control. A 
simple-minded calculation assumes that, for example, the angle ip is changing with time as 



V'(t) = + V'o • t + ^V^o • + 0(t^), (64) 

where tpo and ipQ are the constant rate of this angular change and constant angular acceler- 
ation (a similar analysis could be made for the angle e). Assuming that for any given time 
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the quantity iP{t) has a normal Gaussian distribution with errors obeying the equation: 

AV^(t) = AV'o-t + C(t2). (65) 

the later expression may be taken as an input to derive the on-board attitude requirements. 
Let the time interval for sampling will be of the order of r — At s then from Eqs.(41) and 
(65) one derives the following constraint: 

Ak 

a^^ ~ 38.529 • — mas/s. (66) 

Thus, for the sampling rate of r = 0.5 s (e.q At — 0.5) and Ak — 0.701 (for absolute 
angular acceleration compensation, e.q. -00 = 0), from this last equation (66) one obtains 
the following requirement on the accuracy of spacecraft attitude control: 

crf^l ~ 54.018 mas/s. (67) 

Concluding, we would like to mention that the derived requirements Eqs.(43),(67) were 
based upon the assumption that the total error budget for the velocity aberration is — 
0.701 lias. Naturally, minimizing the error budget for the velocity aberration cTq^ will result 
in tightening the derived constraints on the knowledge of the velocity vector itself v — 
V ( cos -0 COS e, sin cos e, sin e ) and the corresponding ay,a^,ae- 



Conclusions 

We have demonstrated that the 4 mm/s requirement for the relativistic stellar aberration 
is quite well justified for the SIM related studies. However, a secondary effects, such as 
possible correlation between the constant term and the instantaneous spacecraft velocity 
vector may produce a noticeable impact on the astromctric observations. Thus, account for 
this possible correlation, actually results in tightening the velocity requirement (for the worst 
case of highly correlated data) and almost twice minimizes a tolerable error in the velocity 
estimates, notably = 1.616 mm/s. Additionally, the solar radiation pressure will greatly 
influence the SIM orbital solution and will be the major force acting on the spacecraft. As 
a result the navigation DSN time will not go down as the mission progresses because of 
fluctuations in the solar radiation pressure. 

These issues may be addressed by analyzing a number of different options that will 
present a grounds for possible trades between the Mission, Spacecraft, and Instrument Sys- 
tems. Some of these options are 1). On-board processing of data, as oppose to increasing the 
use of DSN time; 2). The use of a precisely positioned solar shade for maintaining constant 
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pressure loading; 3). Velocity determination using integrated and time averaged accelerom- 
eters. In any case the development of a better model for the instrument in order to perform 
analysis of its 'elastic'/dynamical properties and the systematics hidden in the "constant" 
term will significantly boost the related studies. 

Another issue that needed to be discussed is the reducing the velocity accuracy require- 
ment down to 2 mm/s. A first look at this problem brings a conclusion is that this is not 
an easy problem for SIM in deep space (for more details, please see You and Ellis (1998); 
Gorham (1998)). Thus, some of the issues concerning SIM in a SIRTF-like orbit, with a 
2mm/s velocity knowledge requirement are leading to the following conclusions from the 
DSN tracking stand-point: 

1. By utilizing the Doppler range-rate method one can achieve more than adequate 
precision in the radial direction, of order 0.05 mm/s is reasonable. But in the plane of the 
sky, it is more difficult: long observations (doppler arcs) are necessary and it is unclear that 
the same precision can be obtained uniformly in any case - it will likely be a function of the 
ecliptic latitude of SIM at the time of observation. 

2. By using the Delta-DOR observations (VLBI) with a single baseline and X/S-band 
observations, one may be able to obtain a precision of at best ~ 1 nrad for a fairly short 
observation (tens of minutes, including QSO calibrator). This corresponds to ~ 15 m at 
0.1 AU. So, a pair of observations separated by ~ 2 hr or, so could get one plane-of-sky 
component to the required accuracy; ideally a track would be used to improve the SNR 
and get somewhat better accuracy perpendicular to the baseline. To get both plane of sky 
components a pair of baselines is necessary [a bigger DSN commitment, of course]. 

It is unclear whether the radial Doppler range-rate will be available with the same 
antennas at the same time as the VLBI observations. Development of a Ka band capability 
would improve the precision by a factor of ~ 3. This is why meeting the SIM requirements 
is not completely trivial. It might improve if the distance is not as far, and it would be 
interesting to see how performance varied with a number of parameters, including declination. 
There is another issue which was raised in this paper. That is a reasonable definition for the 
data sampling rate At necessary to use in order to introduce an on-board correction for the 
relativistic aberration. Our results show that the answer to this question may impact the 
whole Mission operations approach. This is why we believe that this issue should be a topic 
for a separate study and the corresponding results will be reported elsewhere. 

Let us mention that the resulting SIM astrometric catalog will be derived by reducing 
the collected data to the solar system barycenter. In this case the uncertainty in Earth's 
barycentric velocity may be additional source of errors. Let us estimate the accuracy of the 
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presently known value for Earth' barycentric velocity f®. It is known that the Earth velocity 
uncertainty is dominated by the uncertainty in the Earth's orbital plane orientation, at the 
~1 mas level, giving a velocity uncertainty of ~0.03 mm/s. However, there may be given 
another more conservative estimate to this quantity by using the accuracy of determination 
of the Earth-Moon barycenter, which is presently known to accuracy of about ~ 1 km. This 
number helps to estimate the present accuracy of knowledge the Earth's barycentric velocity 
as = 0.2 mm/s, thus offering a fair possibility for ~ 1 //as astrometry [that, as we saw, 
requires ~ 4mm/s]. So this is not a problem for SIM. 

The requirements on the navigational accuracy obtained here, relate to the interferom- 
eter itself and not to the center of mass of the spacecraft. In the case of changing of the 
center mass of the spacecraft (which may occur due to different reasons, such as the rotation 
during imaging mode, or due to the use of propellant, etc.), the navigational parameters 
should be re-calculated to the required accuracy. Let us also mention that at the SIM level 
of accuracy a more subtle effects on the motion of the Earth about the sun will start to 
play an important role. These effects arc the Moon's and the largest planets influence, and, 
probably, the body-body interaction between the Earth and the Moon, the spin-orbital con- 
nection in the Earth dynamics, the solar rotations and pulsations, influence of solar plasma, 
etc. All of these effects will have to be calculated and properly included into the models for 
astrometric grid simulations and corresponding data analysis. Our further studies will be 
aimed on improvement of the models used for the presented analysis by developing a better 
model for the instrument behavior and by expanding the parameter space to incorporate the 
other physical phenomena affecting astrometric observations from within the solar system. 

The reported research has been done at the Jet Propulsion Laboratory, California Insti- 
tute of Technology, which is under contract to the National Aeronautic and Space Adminis- 
tration. 
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A. General expression for the errors propagation 

In this Appendix we will present the intermediate calculations omitted for brevity in 
the main text. Parameterization Eq.(8) 

si — (cosaicos^i, sincKiCOS^i, sin^i), 

S2 — (cosq;2COS(52, sinQ;2C0s52, sin52), 

V — V ( cos ip cos e, sin ip cos e, sin e ) , 

b — 6(cos/xcosi/, sin cos I/, sini/). (Al) 

allows one to express the scalar products in Eq.(7) in the usual way: 

{b ■ Si) — 6^cosi/cos(5i cos(q;i — //) + sin i/sin , 
{b ■ S2) = 6^ cosi^cos52 cos(q;2 — A*) + sin i/sin (^2^ , 
{v ■ si) = cos e cos Si cos(q;i — ip) + sin e 

5 S2 cos(q;2 — ■0) + sin e sin 5; 



5£ = b 



{v ■ si) = v ^cos e COS 5i cos(q;i — t/^) + sin esin , 
{v ■ S2) — cos e cos S2 cos(q;2 — ip) + sin e sin ^2^ , 
. present equation (7) in the fully-blown form as follows: 



(A2) 



As a result, we can 

— /i)^ + sin 1/ ^ sin 5i — sin 62^ 



cosui 



' ^ COS 5i cos(q;i — jj) — cos 82 cos(q;2 
cos V COS 5i cos(q;i — /x) + sin v sin ^ cos e cos 61 cos(q;i — -0) + sin e sin 81^ — 
— ^ cos V COS 82 cos(a2 ~ A*) + sin v sin 82^ ^ cos e cos 82 cos(q;2 — V') + sin e sin 82^ ■ (A3) 



bv 
c 



By taking the first derivative from both left and right sides of this equation (A8), we 
will obtain expression necessary to analyze the contributions of the different error factors 
involved in the problem under consideration to the overall error budget: 



1=1,2 



fd8t ^ d8e^^\ d8i 

Arv.- -\ AS: 1 -I 



dSe , d8i , d8i , d8i ^ d8i ^ , , , 
+ —A// + -^Au + —Av + —Ae + -^Ai^ + A8cq. 
o/i ov ov oe dtp 



(A4) 



with the partial derivatives with respect to the observing angles 5i, 82 and ai, a2 given as 
follows: 
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m 
882 



dse 



dSi 



— ^ I ~ cos u sin di cos(q;i — /x) + sin u cos Si + 

+ - ^ cos e cos 5i cos(q;i — ip) + sin e sin ^ cos v sin 5i cos(q;i — 11) — sin v cos + 
+ ^ cos e sin 5i cos(q!i — ip) — sin e cos Si^ ^ cos i/ cos di cos(q!i — /j,) + sin 1/ sin di^ j | , 
= 6 I cos v sin ^2 cos(q;2 — jj) — sin 1/ cos 82 — 

— — 1^ ^ cose cos (^2 cos (q;2 ~ V^) + sine sin (^2^ ^cos 
+ ^ cos e sin ^2 cos(q;2 — ip) — sin e cos ^2 j ^ cos cos 
h J — cos cos 5i sin(Q;i — //) + 



5 sin (52 cos(q;2 ~ A*) ~ sin i/ cos 
5 (^2 cos(q;2 — //) + sin u sin 02 



= ^ { ~ 

+ - cose cos (5i cos(q;i 



62) + 
sin52)]}, 



— tp) + sin e sin 5i j cos u cos (5i sin(ai — fj,) + 
+ ^ cos cos 5i cos(q;i — /x) + sin sin cose cos (5i sin(Q;i — -0) |, 



= 6 I cos v cos 62 sin(a2 — 1^) 



— - 1^ ^ cos e cos (^2 cos(q;2 — "0) + sine sin ^2^ cos cos (^2 sin(Q;2 — A*) + 
+ ^ cos v COS (^2 cos(q;2 — //) + sin 1/ sin ^2^ cos e cos ^2 sin(Q;2 — i^)^ | , 



(A5) 



together with the partial derivatives with respect to components of the basehne vector b: 
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— cosz/( 



+ 



^ COS 6i cos(q!i — jji) — COS ^2 cos(q;2 ~ y"-)^ + sin v ^ sin 5i — sin + 
OS e cos 82 cos(q;2 — ■0) + sin e sin ^2 j ^ cos v cos ^2 cos(q;2 — /i) + sin sin ^2^ — 

— ^ cos e cos 5\ cos(q;i — ip) + sin e sin 5i j ^ cos v cos 5i cos(q;i — /x) + sin 1/ sin , 
= 6 I sin ^ cos ^2 cos(q;2 — jj) — cos 5i cos(q;i — ji)^ + cos ^ sin 5i — sin ^2^ + 

— ^ cos e cos ^2 cos(q;2 — 'V') + sin e sin ^2 j ^ — sin i/ cos 82 cos(q;2 — A*) + cos v sin 82^ — 

— ^cos ecos^i cos(q;i — ■0) + sinesin^i j ^ — sin i/cos5i cos(q;i — jj) -\- cos sin 5i 



V 

+ - 
c 



d5l 



= b cos u 



V 

+ - 

c 



'I COS 61 cos(q;i — jJi) — cos 62 cos(q;2 — A*) + 

s e cos 82 cos(q;2 — V') + sin e sin ^2^ cos 82 sin(Q;2 — jj) — 
— ^ cose cos 5i cos(q;i — 'i/^) + sine sin 5ij cos 5isin(Q;i — /i) |, 



(A6) 



and, finally, with the following partials for the spacecraft's velocity: 



dv 



d8i 



cos e cos 82 cos(q;2 — ip) + sin e sin 82^ ^ cos u cos 82 cos(q;2 ~ A*) + sin 



1/ 



sin 82 



— ^ cos e cos 81 cos(q;i — ■0) + sin e sin ^1 j ^ cos 1/ cos 5i cos(q;i — //) + sin 1/ sin 



— — = — cos e 



in "ly ^ cos // cos ui cos\^U!i — /j,) -t sin u sin oi^ , 

sine COS (5i cos(q;i — ip) — cos e sin 81^ ^ cos 1/ cos 81 cos(q;i — //) + sin 1/ sin — 
82 cos(q;2 — i^) — cos e sin ^2 j ^ cos u cos ^2 cos(a2 — /i.) + sin f sin ^2 j , 
cos V COS (^2 cos(q;2 — //) + sin 1/ sin ^2^ cos 82 sin(Q;2 — ip) — 
^ cos i/cos^i cos(q;i — /x) + sin i/sin 5i j cos 5i sin(Q;i — ■0) . (A7) 



bv 
c 

— ^ sin e cos 82 
bv 



The group of the expressions Eqs.(A4)-(A7) is quite difficult for analytical description, 
however it may be significantly simplified. First, without loosing generality one can set 
Ij, = 1/ = 0, which equivalent to choosing the direction of the baseline vector b coinciding 

— * 

with X-axis, namely 6 = 6(1, 0, 0). As a result of this choice, all vectors now will be counted 



-33- 



51 



hv 



COS( 



from the baseline and the expression (A3) for the relative OPD may be presented as follows: 

cos 5i cos Oil — cos 82 cos + Scq — 

5 5i cos tti ^ cos e cos 5i cos(q;i — ip) + sin e sin — 

— cos 82 cos q;2 ^ cos e cos ^2 cos(q;2 — ■0) + sin e sin ^2^ . (A8) 

Second, the largest expected ratio ^ w ^ = 9.942 x 10~^, which makes the terms ~ ^ in the 
expressions (A5) and (A6) of the second order of smallness. For the purposes of the present 
analysis, we may omit these terms and, after some re-arranging, expression Eq.(A4) may be 
presented as 

— — — Aq;2 sm a2 cos 02 — Acti sm cti cos Oi + 




+ A(52 cos a2 sin 82 — A5i cos cci sin 5i + 

, ^r...., .... ..... ....1 , ^'^^o 

Av 



cos QJi COS 5i — COS q;2 cos 82^ 



+ —r^ + 



+ — 
c 



COS e ( cos^ 82 cos q;2 cos(q;2 — ip) — cos^ 5i cos a\ cos(q;i 



+ 



+ Ae - 



+ sine ^ cos (^2 sin (^2 cos q;2 — cos 5i sin 5i cos + 

sin e ^ cos'^ (5i cos ai cos(q;i — ip) — cos^ 82 cos 0:2 cos(q;2 — + 

(^2 sin (^2 cosq;2 — cos^i sin (5i coscci^ + 

cos a2 cos^ ^2 sin(Q;2 — — cos cti cos^ 5i sin(Q;i — tp) . 



+ cos e ^ cos (^2 ! 



+ A'^ - cos e 
c 



(A9) 



This last expression may further be simplified by introducing a parameterization, which 
is natural for the problem under consideration. The final result for the first derivative of the 
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differential OPD may be given as follows 
A5£ 



^ Ace sin a-2 cos 82 + Acti ^ sin q;2 cos 82 — sin ai cos + 

+ A5 cos 0:2 sin ^2 + A5i ^ cos ^2 sin 82 — cos cti sin 61^ + 
A6 / ^ \ A5co 

+ — cos «! COS Oi — COS a2 COS ()2 j H ^ h 

'"^^ cos(e — eo) — Ae - sin(e — eo)j ^J {a? + p) cos^('0 — V'o) + 



+ 



— A'0 - cose-\/a2 + Psm^ij) — i/jq). 



(AlO) 



where both angles ipo and eo are depend only on the positions of the two stars involved and 
are given as follows: 



sin ipQ 

a 
f 

sin eo 

cos eo 
k 



f 



cos ijjo 



cos^ 82 cos^ a2 — cos^ 81 cos^ ai, 

cos^ 82 cos a2 sin «2 — cos^ 81 cos cti sin ai . 
k 



(All) 



v/(a2 + /2)cos2(V'-Vo) + A;2' 

cos 82 sin (^2 cos 0:2 — cos 81 sin 5i cos ai . 



(A12) 
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B. Errors in the fully-parameterized relative OPD 

Taking into account the fully parameterized form of the fractional relative OPD Eqs.(A4)- 
(A7) we may obtain the form of this quantity for a tile in the SIM nominal observing direction. 
Thus, substituting in these expressions the values for the positions of primary and secondary 
stars Eq.(14) we will have: 



— r— = Aa cosi/cos - cos(— — /x) + 2A(q;i — ^) cosi/cos - sm — sm// 



- AS 

- 2Aiy 
2Ab 



cosi/sm-sm( //) + smi/cos- + 2 AOi cos sm - cos — sm // 

2 2 Zi - 2 2 



6 . a 



6 



+ 
+ 



b 

Av 
c 



sm v cos - sm — cos a + cos v sm 

2 2^ 2J 

S a S 
cos u cos - sm — cos a — smu sm - 
2 2^ 2J 



+ 



ASco 
+ —r^ + 



a 



cos V I — cos e sin a cos - sin(^ + /x) + sin e cos — sin 6 sin fi 



+ Ae 



+ sin u cos e cos ^ sin S sin ■^j + 
cos u(^sme sin a cos^ - sin('0 + /x) + cos e cos ^ sin 5 sin //^ — 



q; 



— sin u sin e cos — sin 5 sin -0 
2 ^. 



+ 



+ A-^ -cose ( sini/ cos — sin5cos'0 — cosi/sincKCOS^ - cosfV' + ). (Bl) 
c v 2 2 / 



One may see that the angles of the baseline orientation {fi, v) are significantly infiuencing 
the narrow-angle astrometric observations. However, as previously we will choose both angles 
as /i. = = 0, which is equivalent to choosing the direction of the baseline vector b to coincide 

— * 

with X-axis, namely b — 6(1, 0, 0). As a result of this choice, all vectors now will be counted 
from the baseline (see Figure 2). Resulted expression for the relative OPD may be presented 
in a simpler form, namely: 



A51 , a 5 
—— — Aa cos — cos - 
b 2 2 

'Av 



AC- " ■ ^ 

Ad sm — sm - 
2 2 



2A6 a 5 A5co 

sm — cos - + — ; — 

b 2 2b 



( — - cose — Ae - sine) sin'^ + A'^ - cose cos V'] sincccos^ -. (B2) 
\ c c / c J 2 
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Note that obtained expression does not depend on the errors in position of the primary star 
Aai, ASi. The obtained result may now be used to analyze the propagation of errors in the 
future astrometric observations with SIM. 
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c. Components of the covariance matrix 

The first term in the expression (B2) A(5£ may equivalently be presented as A5f = 
A(nAo) = An Ao + n AAq, where Aq is the operating frequency and n is an integer number. 
This term vanishes because both Aq and n are assumed to be known to a sufficiently high 
accuracy, e.q. An = 0, AAq = 0. This is true due to the fact that the SIM white hght 
fringe tracker will make sure that variations in these two quantities will be exactly zero. 
Then the remaining differentials Ai>, Ae, A'0, Aa, A5, A6 and A5co will have to satisfy 
the equation: 



A a 5 
Aq; cos — cos - 
2 2 



... « . 5 A5co 

Ao sm — sm - = ; — 

2 2 h 



2Ab . a 6 
— — sm — cos — h 
6 2 2 



+ 



/Av 
— cose 
\ c 



— Ae - sin e I sin + A^^ - cos e cos ijj 
c / c 



sm a cos 



(CI) 



In order to simplify further analysis we have separated the terms in the expression above in 
a such a way, so that the left side of this equation represents the error in the measurement 
of the absolute angular separation between the two stars, while the left side shows the main 
contributing factors to this quantity. 

In order to study the error propagation in the astrometric observations, we need to take 
the square of the expression (CI): 

/A\2 2*^ , ^ r^2 ■ 2 ■ 2 ^ (AaA6) . . . 

(Aa) cos —cos - + Ad sm — sm smctsmo = 

^ ^ 2 2^^ 22 2 



+ 



jASco] 

Av 
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2ASco 
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4Ab 



+ 



cose 



4(A6)2 . 2^ AA5coAb) . a S 

sm^ - cos - + 4^ sm — cos - + 

62 2 2 62 2 2 



— Ae - sin e I sin + A'^ - cos e cos ijj 
c J c 



2 S 
• 2 4 " 

sm a cos - 
2 



/Av 
— cose 
\ c 



A 

Ae - sm 
c 



in e^ sin ip 



V 



+ Alp - cose cos 
c 



/Av 
\ c 



cose 



Ae - sin e ) sin tjj + Aijj - cos e cos tjj 
c / c 



sm a cos 



a 



smasm — cos -. 
2 2 



(C2) 



One may expect that in any given tile the errors in Aa and AS are normally distributed 
and uncorrelated.^ The errors due to the orbital dynamics Av, Ae, A^^ at the chosen 



^Note that this is not true for a general case of studying the reference grid stabiUty. Thus one finds that 
the correlation in RA and DEC becomes a source for the zonal errors in the analysis of the grid accuracy. 
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approximation may also be treated as being not correlated with instrumental errors Ab, Acq. 
As a result, the components of the covariance matrix may be given as follows: 



{Aay 




{ASf 




AaA5 


= 0, 


{Avf 




{Aef 








AvAe 


= 0, 


AvAip 


= 0, 


AeAip 


= 0, 


{AScoy 




A5c,Av 


= 0, 


(A6)2 




AbAv 


= 0, 


AbAe 


= 0, 


AbAi/j 


= 0. 



(C3) 



A6coAb = (JScn<^bP{cQ,b), A6coAe = as o-^p{co,e), A6coAij = as p{cq,^). 



(C4) 



These expressions are helpful to present the result of averaging the equation (C2) in the 
following form: 



(T„ COS — COS --\-ai sm 
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sm 



62 
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4 a, 



^ sm - cos - + — p — p(co, 6) sm 2 cos - + 



+ 



+ 



a. 
6 



2 2^ •2^*2/ 2^ 2 2/1*2 4*^ 

COS e + (j^ ~2 ^ j W ^ ~^ cos e cos -i/^ sm a cos - + 
(jgp(co, e) sinesin'^ — cr,^ p(co, "0) cos e cos - sin a cos^ -. (C5) 



Expression (C5) suggests that for the observations in the direction perpendicular to the 
baseline (e.q. when 5 are non-zero), there will be a large contribution to the error budget 
coming from as- However, for observations parallel to the baseline (e.q S — 0) the influence 
of as is vanishes and one obtains the most stringent requirement on the velocity aberration. 
Then, by taking S — one may present Eq.(C5) in the following form: 
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/ C"^ 2 2^ •2^*2/ 2^ 2 2/ 
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(Te p(co, e) sine sin — p{co,'ij:) cose cos iIj 



V 

- sma. 
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(C6) 



Remembering further that angles {a, 6) vary in the range given by + 5^ < (f^)^ and cos ^ 
never vanishes in this interval, we can divide the both sides of the equation (C6) on cos^ |: 



62 cos2 I 



+ 



62 2 + 



62 



COS 



2 a 



+ 4 

4(T5, 
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f(^v 2 . 2 V . 2 \ ■ 2 I . 2 V 2 2/1.2", 
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- tan 
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(C7) 



This equation represents an ellipsoid of with half-axis defined by the direction of motion 
of the spacecraft. Thus, for the spacecraft motion in the direction given by = f and e = 
we can have the first expression for astrometric errors: 



. _ <o , 4^2 a ,^ sin f 4(72 a 

- 6^W^ + "6^ 2 + —6^^^"°' ^ + ^ 2 ■ ^^^^ 

For the case of motion in the direction = f and e = ±| we have the second relation: 



2 ^lo , 4a2 a 4a, sm f 

6Wf + -6^'"" 2 + -6^^^'°'^^^ + 



+ 4a2 — sm2- ± _ p(co, e) tan-. (C9) 



And, finally, the last expression may be obtained when considering the motion in the direction 
■0 = and e — 0, tt. This last constraint reads: 

2 ^5cn 4a? o q; 4a<5 ab sinf 

+ 4a; — sm^- ^ __JLa^_p(co,V') tan-. (CIO) 
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For the simplicity, we can relate the errors in the sky angles of velocity Ue and with that 
of the velocity magnitude cr^. To obtain the most stringent constraints on the sky angle we 
will use to following expressions: 



a 



^^p(co,e) tan- 



V 



a 



a. 



a 



^ p(co,^)tan- 



(Cll) 



As a result, in our further analysis we will concentrate on the equation Eq.(C8) together 
with the following two solutions for and (which were obtained directly from Eqs.(Cll)): 



, r^^co p(co,e) c]2 (Ts,^ p(co,e) c 



+ 



sm O' ? ' 



sma V 



> 0, 



sm a i> 



> 0. 



sm a V 



(C12) 



Just for the estimation purposes we can assume that p{co,ip) = p(co,e) = po- Then, by 
using the expressions (C12) together with the estimates Eqs.(35) one could obtain improved 
requirements on the accuracy of knowledge of the two sky-angles of the spacecraft velocity: 



(7e = (7v, = a«--^(yi+^-po) >0. (C13) 
V 2 sm f V V / 

Note, that for the worst case of highly correlated quantities e) and Cq, e.q. po ~ 1 the 
result Eq.(43) will have to be further reduced by a factor ^/2 — l — 0.4142, thus tightening the 
requirements for the accuracy of knowledge of the two velocity sky-angles as a^ — a^ — 20.186 
mas. This example suggests that a possible correlation between the velocity components and 
the constant term may be reduced if the DSN will deliver the data for all three components 
with a twice higher accuracy then planned. 



